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Abstract 

Thermodynamically stable vortex-antivortex structures in a quasi-two-dimensional superconduc- 
tor in a tilted magnetic field are predicted. For this geometry, both orbital and spin pair-breaking 
effects exist, with their relative strength depending on the tilt angle 9. The spectrum of possi- 
ble states contains as limits the ordinary vortex state (for large 9) and the Fulde-Ferrell-Larkin- 
Ovchinnikov state (for 9 = 0). The quasiclassical equations are solved near H C 2 for arbitrary 9 and 
it is shown that stable states with coexisting vortices and antivortices exist in a small interval close 
to 9 = 0. The results are compared with recent predictions of antivortices in mesoscopic samples. 
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I. INTRODUCTION 



In this paper we shall establish a link between two current topics in superconductivity 
which seem unrelated at first sight. These two topics are antivortices in mesoscopic samples 
and paramagnetic pair breaking effects in layered superconductors respectively. 

Antivortices are topological singularities of the order parameter with vorticity oppo- 
site to ordinary vortices. Generally, the phase x °f the superconducting order parameter 
ip = \ip\expzx must change by 27m, when surrounding an isolated zero of ip which is a 
topologically singular point. The arbitrary integer n is referred to as topological charge 
or vorticity. An ordinary vortex in a type II superconductor has positive vorticity, n — 1, 
corresponding to a phase change by 2n, a total flux of a single flux-quantum $o, an d a dia- 
magnetic behavior of the (screening) currents around the vortex center. An antivortex, on 
the other hand, has negative vorticity (in the simplest case n — — 1) and carries a negative 
flux (e.g. — <&o)- The currents surrounding the antivortex center have opposite direction as 
compared to the ordinary vortex, i. e. they are paramagnetic in nature. The diamagnetic 
nature of the screening currents in the ordinary, orbitally-limited, superconducting state is 
a consequence of the dominating orbital pair-breaking mechanism. 

Notwithstanding this fundamentally unfavorable conditions for the occurrence of an- 
tivortices, several specific configurations designed to create stable vortex-antivortex struc- 
tures in orbitally-limited superconductors, either in thermodynamic equilibrium 1 ' 2 ' 3,4 or as 
metastable states 5 , have been proposed in recent years. Of particular interest are theoretical 
predictions in mesoscopic superconducting samples 1,A . Free energy calculations in the frame- 
work of GL-theory, following the original prediction 1 of stable vortex-antivertex structures, 
revealed 6 ' 7 that such structures are in fact unstable in type II superconductors; this result 
was obtained using a GL-parameter k considerably larger than the critical value k = 1/ \[2 
separating type II from type I superconductivity. This instability is a consequence of the 
attraction between vortices and antivortices which consequently annihilate each other in the 
state of lowest free energy. It may be qualitatively understood by considering the interaction 
between well separated vortices in bulk samples 8 . In this case, the total interaction between 
vortices in GL theory is the sum of two contributions, a repulsive electro dynamic part and 
an attractive condensation energy part. The electrodynamic part exceeds the attractive 
part for k > 1/ \/2. But for a vortex-antivortex pair the electrodynamic interaction changes 
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sign and becomes attractive while the condensation energy interaction does not change sign 
and remains attractive. Thus, the total force between vortices and antivortices becomes 
attractive. 

Nevertheless, the search for stable vortex-antivortex structures in mesoscopic samples 
was finally successful 4 . A state with total flux of 2$ , carried by a antivortex and three 
ordinary vortices in a mesoscopic sample of triangular shape, was found to be stable in a 
small region below T c and below k = l/y/2 (In this part of the «, T— plane, GL-theory still 
applies 9 ). This effect, which has not yet been verified experimentally, is essentially a surface 
effect and cannot be understood in terms of forces between isolated vortices. Contrary to the 
interpretation of the authors 4 , the decrease of k below the critical value l/y/2 does not lead 
to a change in sign of the vortex-antivortex interaction. The latter remains attractive for 
arbitrary k since only the electrodynamic part of the interaction changes sign. Nevertheless, 
the decrease of k plays an essential role because it weakens the relative importance of 
the vortex interactions in comparison to that of the confining boundary conditions. The 
geometric shape (triangular) of the sample is important since no stable vortex-antivortex 
structure has been found 10 in thin films of quadratic shape. 



II. PARAMAGNETIC PAIR-BREAKING AS A STABILIZING FACTOR FOR 
ANTIVORTICES 

Considering the physical mechanism governing the interaction between vortices it seems 
impossible to create periodic structures with antivortices for values of k much larger than 
1 / \/2, deep in the type II region. This is in fact true for the purely orbitally-limited supercon- 
ducting state. However, as will be shown in this paper, an exotic type of superconductivity 
exists, characterized by the presence of two different competing pair-breaking effects, where 
vortex-antivortex structures appear quite naturally for a large range of GL parameters k, 
not as a surface effect but in an extended, two-dimensional (2D) periodic "vortex" -lattice. 

Cooper pairs may either be broken by the orbital effect or by means of the interaction 
between the external field and the magnetic moment due to the spins of the electrons 11,12 . 
This spin pair-breaking effect leads to a paramagnetic response of the superconducting con- 
densate. It is usually much smaller than the orbital effect but becomes important if the 
orbital effect can be suppressed. Clearly, decreasing the magnitude of the orbital effect 



3 



favors the stability of antivortices. 

Until recently, in most works only the limit of complete suppression of the orbital effect 
was considered. In this paramagnetic limit, a spatially inhomogeneous superconducting 
state, referred to as FFLO state, has been predicted theoretically by Fulde and Ferrell 13 , 
and by Larkin and Ovchinnikov 14 . With regard to the present problem, the FFLO state 
itself is not a good candidate for observing antivortices since the diamagnetic response is 
completely suppressed, and no (anti)vortices at all can exist in the purely paramagnetic 
limit. 

Obviously, coexistence of both pair-breaking mechanisms, with only partial suppression 
of the diamagnetic effect, provides optimal conditions for the stability of vortex-antivortex 
structures. Such a situation may be realized by slightly tilting the applied magnetic field, by 
an angle 9, from the plane-parallel direction . If the conducting layers are completely decou- 
pled and infinitely thin (it has been shown recently that single-atomic layers are required for 
negligible orbital pair-breaking 15 ) one has, in fact, a mixture of both pair-breaking mech- 
anisms, with the orbital effect entirely due to the perpendicular component H± = Hsin9 
and the spin effect mainly due to the parallel component H». The relative weight of the two 
pair-breaking effects may be controlled with the help of the angle 9. 

The second order phase boundary H c2 for this situation has been calculated by Bu- 
laevskii 16 for T = and by Shimahara and Rainer 17 for arbitrary T. They found a non- 
monotonic upper critical field curve with different pieces of the curve belonging to different 
Landau quantum numbers n (n — 0, 1, . . .). For large 9 one recovers the usual orbital result 
n — 0. With decreasing 9 unusual pairing states with nonzero n — 1,2, .. . appear for 9 
smaller than a critical angle of the order of one degree. For 9 — > 0, one obtains n — > oo and 
the upper critical field of the FFLO state is recovered. 

III. OUTLINE OF CALCULATION 

The states below the non-monotonic upper critical field curve which belong to the lowest 
nonzero quantum numbers, say n = 1 — 4, are most interesting from the present point of 
view of antivortices. 

We shall make only two very general (and reasonable) assumptions with regard to the 
structure of these states. The first is, that the order parameter modulus |A| and the local 
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magnetic induction B are periodic with respect to an arbitrary two-dimensional lattice. 
The second assumption is that a unit cell of this quasi-periodic lattice carries a single flux 
quantum $ . To find the equilibrium structure of these states, all possible lattices must be 
considered and the global minimum of the free energy must be identified. 

Considering the vicinity of the upper critical field, where a second order transition occurs, 
a free energy expansion, up to terms of fourth order, may be performed. Then, the general 
scheme of such a calculation is the same as Abrikosov's classical work on type II supercon- 
ductivity 18 . However, microscopic equations should be used here, because a GL-formulation, 
with a finite number of gradient terms, is not valid in the present (low-temperature) situ- 
ation. The quasiclassical theory of superconductivity, generalized with regard to Zeeman 
coupling terms, provides an appropriate theoretical framework 19 for the present problem. 

The nonlinear transport equations for a clean superconductor, including Zeeman coupling, 
are given by 

2u>i + hv F (k)'& r f(f, k, Qi) = 2A(r)g(f, k, u t ), 

2Q l -hv F (kK\ f + (r,k,cu l ) = 2A*(r)g(r,k,Cu l ), 

where /, / + , g are the quasiclassical Green's functions (g 2 = 1 — ff + ), A is the order pa- 
rameter, and i? r is the gauge-invariant derivative, defined by $ r = V r — i(2e/hc)A. The 

— # 

latter term, containing the vector potential A, describes the familiar orbital coupling be- 
tween Cooper pairs and magnetic field. In addition the coupling between the spins of the 
electron and the magnetic field is taken into account in Eq. (1). The corresponding Zeeman 
term fiB is contained in the complex variable Qi — out + %[iB which replaces the real Mat- 
subara frequency uj\. The quantity v F {k) denotes the Fermi velocity which depends on the 
quasiparticle wave-number k. 

The order parameter A is defined in terms of the Green's functions /, / + by the relation 

2nk B T J2~ + ln ( T / Tc ) A ( f ) = 7r ^ T E f d2k ' [fP, *i) + f(r, *!)] , (2) 
i=o UJl ) 1=0 ^ 

where Nd is the cutoff index for the Matsubara sums. To calculate the vector potential A 
for given Green's function g, Maxwell's equation 

V r x (B(r)+AnM(r)) = * ^J^f "4 — v F (k )Qg(r, k , u>i), (3) 

c 1=0 71 
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has to be solved. Here, N F is the normal-state density of states at the Fermi level. The 
r.h.s. of Eq. (3) is the familiar (orbital) London screening current while the magnetization 
M is a consequence of the magnetic moments of the electrons and is given by 



The first term on the r.h.s. of Eq. (4) is the normal state spin polarization. The second 
term is is a spin polarization due to quasiparticles in the superconducting state. 

Solving these equations one obtains extrema of a free energy functional G (not written 
down here, see 19 ), whose global minima represent the stable states. The quasiclassical 
equations (l)-(4) are rather general; they describe both orbital and paramagnetic pair- 
breaking effects, and cover the ordinary vortex state as well as the FFLO state in the 
appropriate limits. In addition, Eqs. (l)-(4) take into account the magnetic response of 
the superconductor (arbitrary values of the GL-parameter k may be considered), which is 
neglected in most treatments of paramagnetic pair breaking. Note that the familiar London 
screening current is completely absent in the purely paramagnetic (FFLO) limit. 

The calculation generalizes techniques used previously to calculate the upper critical 
field 16 , 17 and the equilibrium structure in the high-K-limit 20 . For the vector potential the 
same gauge as in previous numerical calculations 9,21 on the vortex lattice (without Zeeman 
coupling) may be used. A first important step is the solution of the transport equation (1) 
for small |A|, taking derivatives of |A| of arbitrary order into account . This is achieved 
by means of an eigenfunction expansion and the so-called "Helfand-Werthammer Trick", 
leading to an integral representation for the first-order Greens functions. In the course of 
the following free energy expansion, up to terms of fourth order in |A|, a large number of 
momentum- and configuration space integrations have to be performed. The evaluation of 
these integrals may be greatly simplified by introducing the gap correlation function. A 
more detailed description of the lengthy calculation (and a discussion of other results, not 
related to antivortices) may be found elsewhere 19 . 

IV. STABLE VORTEX- ANTIVORTEX STRUCTURES FOR N=2 

In this section results are reported which show that stable vortex-antivortex structures do 
in fact exist in bulk superconductors for nonzero finite n. The most important of these states, 
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which is most stable and most easily accessible experimentally, belongs to Landau quantum 
number n = 2. Discussion will be restricted here to this single state (Stable structures for 
other n and a discussion of the general properties of the new paramagnetic vortex states have 
been reported elsewhere 19 ). The results depend on 4 input parameters, the dimensionless 




FIG. 1: Square of modulus of order parameter \ip2 1 2 for Landau level n = 2 as a function of 
x/a, y/a in the range 0.25 < x/a < 1.25, —0.15 < y/a < 0.7. This is the stable structure (unit 
cell parameters a/b = 1, a = 60 ° ) for fj, = 0.1, R = 10, = 0.7 ° , t = 0.2. 

magnetic moment \x = nkBT c /mvp, Eilenberger's GL-parameter k (which is defined in terms 
of the GL-parameter k of a clean superconductor by the relation k = 0.6837k), the tilt angle 
6, and the dimensionless temperature t = T/T c . Given these numbers, the critical field 
B c2 , Landau's quantum number n, the unit-cell structure of the stable lattice (which is 
characterized by the lengths a, b of the unit-cell vectors and by the angle a between them), 
and the fields ip n {r)i B{f) may be calculated 19 . The square of the order parameter modulus 
for pairing in Landau level n is given by the Fourier series 

m 2 (r) = J>D^< (5) 
id 
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with reciprocal lattice vectors Qij corresponding to the lattice parameters a, b, a, and with 
Fourier coefficients 

= (-l) lj e-*< c ° sa e-^/ 2 L n (x ltj ). (6) 

The quantity xij is defined elsewhere 19 . The magnetic response of the superconductor is 
characterized by means of the deviation B\ — B — B of the local induction from the macro- 
scopic induction B. The parallel component of Bi, which describes a purely paramagnetic 
behavior, is denoted by Bn. The perpendicular component B\± describes a mixed orbital 
and paramagnetic response. The Fourier coefficients of these quantities are given by 

t(|A ra | 2 ) 

{Bl\\)l, m ~ ~ ~ 2 _ ^2 B \W ' 



l±)Lm 



-ly^-^^Mx^s^, (?) 

k 2 — fl 2 

(B^ 2 h{xi,m)S\% - giMS^) , (8) 

where (| A n | 2 ) denotes the spatial average of the square of the order parameter and the func- 
tions fx, g% and the Matsubara sums Sj^ are reported elsewhere 19 . Choosing parameters 





FIG. 2: Spatially varying part of the parallel magnetic field (f) using the same input parameters 
and the same part of the x, y-plane as in Figure 1. 



S 



FIG. 3: Spatially varying part of the perpendicular component of the magnetic field Bi±(r) using 
the same input parameters and the same part of the x, y-plane as in Figure 1. 

fx = 0.1, k = 10, 8 = 0.7° , and t = 0.2 one obtains a Landau quantum number n = 2 
(B C 2 = 4.27086 if measured in Eilenberger units 19 ) and a unit cell with a/b — 1, a — 60 ° , 
i.e. a triangular lattice. The order parameter has three zeros per unit cell as shown in 
Figure 1. Thus, one expects that this stable structure contains two vortices and a single 
antivortex, because the total flux per unit-cell is $o- 

The nature of these order parameter zeros may be further clarified with the help of the 
local magnetic field. We plot the fields B\\\ and Bi± as given by equations. (7)-(8), omitting 
a common factor t(|A n | 2 )/ (k 2 — /i 2 ) and using the same input parameters and the same part 
of the x, y-plane as in Figure 1. Let us discuss first the parallel component Bu. As shown 
in figure 2, the functional form of Bu is nearly identical with that of ip2 plotted in Figure 1. 
This behavior, though unexpected at first sight, is a consequence of the purely paramagnetic 
pair-breaking mechanism for the parallel component. The latter implies an enhancement of 
the local induction at points of enhanced order parameter, which is opposite to the familiar 
diamagnetic field expulsion. 

The spatial dependence of the perpendicular component B 1± , plotted in Figure 3, shows 
that the order parameter zero on the left hand side of Figure 1 is an antivortex while the 
two remaining zeros belong to ordinary vortices. At the centers of the ordinary vortices the 
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field has maxima, i.e. shows the familiar behavior of orbitally-limited superconductors. At 
the center of the antivortex the field has a minimum, i.e. shows the opposite, paramagnetic 
behavior. Thus, stable vortex-antivortex structures may indeed be produced by means of a 
2D-superconductor in a tilted magnetic field. It is the simultaneous presence of both pair- 
breaking mechanisms which allows the coexistence of vortices and antivortices in a periodic 
lattice. Calculations 19 for k = 0.1, 1, 100 and t = 0.2, 0.5 show that the structure displayed 
in Figures 1-3 is stable for a large range of parameters. 

Comparing the present results with previous predictions, one notes that the present unit 
cell has the same shape as the mesoscopic sample used in the work of Misko et al. 4 ; the 
number and type of order parameter zeros is, however, different. A second common fea- 
ture is the strong interaction (small distance) between anti/vortices. The basic mechanism 
responsible for the stability of antivortices is, however, completely different in both cases. 
Experimental verification of the structure predicted here requires similar conditions as for the 
FFLO state 22 ' 23 ' 24 ' 25 , namely a layered superconducting material of high purity with nearly 
decoupled conducting planes and a very accurate adjustment of the direction of the applied 
magnetic field. Appropriate materials include the intercalated transition metal dichalco- 
genide TaS 2 — (pyridine) , the organic superconductor k — (BEDT — TTF) 2 Cu(NCS) 2 , or 
the magnetic field induced superconductor A — (BETS) 2 FeCl4. The predicted phenomena 
do not sensitively depend on k and /i but large values of these parameters are favorable, 
because of an associated enlargement of the paramagnetic region. A non-monotonic phase 
boundary in a tilted magnetic field has been observed recently 25 in the heavy-fermion su- 
perconductor CeCoIn^. These data may be related to the present prediction, despite the 
fact that CeCoIn 5 differs considerably from the simple model considered here. 
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